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ABSTRACT 

Many bacterial genes are controlled by metabolite 
sensing motifs known as riboswitches, normally 
located in the 5' un-translated region of their 
mRNAs. Small molecular metabolites bind to the 
aptamer domain of riboswitches with amazing 
specificity, modulating gene regulation in a 
feedback loop as a result of induced conformational 
changes in the expression platform. Here, we 
report the results of molecular dynamics simulation 
studies of the S-adenosylmethionine (SAM)-II 
riboswitch that is involved in regulating translation 
in sulfur metabolic pathways in bacteria. We show 
that the ensemble of conformations of the unbound 
form of the SAM-II riboswitch is a loose pseudoknot 
structure that periodically visits conformations 
similar to the bound form, and the pseudoknot 
structure is only fully formed upon binding the 
metabolite, SAM. The rate of forming contacts in 
the unbound form that are similar to that in the 
bound form is fast. Ligand binding to SAM-II alters 
the curvature and base-pairing of the expression 
platform that could affect the interaction of the 
latter with the ribosome. 

INTRODUCTION 

Antibiotic resistance is viewed as one of the most pressing 
public health problems around the world. According to 
the US Food and Drug Administration, the need for 
new antibiotics is imperative to the successful treatment 
of new resistant bacterial infections (1). Several major 
diseases caused by microbes, such as Tuberculosis, are 
becoming increasingly more difficult to treat with 
commonly used drugs. The excess use of antibiotics and 
evolutionary adaptations of bacteria have exacerbated this 
problem. The need to develop new antimicrobial drugs is 
of the utmost importance, and a greater understanding of 
bacterial gene regulation is a key to solving this problem. 



One possible drug target for the new class of antibiotics 
is the recently discovered metabolite sensing domain of 
bacterial mRNA (2-5). About 4% of all bacterial genes 
are regulated at the post-transcriptional level by small 
metabolites binding to mRNAs (6-9). These transcribed 
mRNAs typically have a ligand-binding 'aptamer' 
domain and an 'expression platform' located in the 5' 
un-translated region (UTR). The binding of cognate 
metabolites, which are generally synthesized by the 
proteins of the respective mRNAs, to the aptamer 
triggers conformational changes in the expression plat- 
forms resulting in initiation or termination of either tran- 
scription or translation. The structured 5' UTR, known as 
a 'riboswitch', effectively functions as an environmentally 
sensitive switch capable of sensing local cellular concen- 
trations of metabolites or other ligands in a feedback loop. 
At the level of translation, riboswitches act by exposing or 
obscuring the Shine-Dalgarno (SD) sequence and the ini- 
tiation codon, which are vital elements in ribosomal 
binding and processing of mRNAs (10). When the bases 
in the expression platform are hidden, the ribosome 
cannot bind to the mRNA or initiate translation (11,12). 
When they are exposed, the ribosome can bind and initiate 
or resume translation. 

Many riboswitches that regulate a host of biosynthetic 
pathways have been discovered (13). The metabolites they 
bind include glycine (14), lysine (15-17), adenine (18,19), 
guanine (20,21), 7-aminomethyl-7-deazaguanine (22,23), 
thiamine pyrophosphate (8,24-27), flavin mononucleotide 
(28,29), vitamin Bi 2 (30,31), glucosamine-6-phosphate 
(32-35), and S-adenosylmethionine (SAM) (36,37). SAM 
is a co-factor involved in the ubiquitous methionine bio- 
synthesis and related sulfur metabolic pathways in all 
bacteria (36). Interestingly, the overall fold of the SAM 
sensing riboswitch is different in different types of 
bacteria. Five different folds of SAM riboswitches have 
been discovered to date, vis-a-vis, SAM-I (36,38), 
SAM-II (39,40), SAM-III (41), SAM-IV (42) and 
SAM-V (43,44). Of these riboswitches, SAM-I and 
SAM-IV share similar binding pockets and significant 
structural elements (42), and SAM-V appears to share 
similar core structure to SAM-II (13). SAM-II riboswitch 



*To whom correspondence should be addressed. Tel: +1 404 413 5564; Fax: + 1 404 413 5505; Email: dhamelberg@gsu.edu 
© The Author(s) 2009. Published by Oxford University Press. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/ 
by-nc/2.5), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. 



Nucleic Acids Research, 2010, Vol. 38, No. 4 1393 



is found mainly in proteobacteria and is very specific for 
SAM (10). SAM-II riboswitch binds SAH (5-adenosylho- 
mocysteine) and other similar byproducts of SAM with 
much lower affinity (40). The specificity of SAM-II for 
SAM and the fact that the X-ray crystal structure of 
SAM-II is one of the only riboswitches containing all of 
the sequence of the ligand-binding domain and the regu- 
latory expression platform make SAM-II an excellent 
model system to study (10). 

Furthermore, the crystal structure of the SAM-II 
metabolite-bound riboswitch, solved by Gilbert et al. 
(10), is a classic H-type pseudoknot shown in Figure 1 
(45), a motif that is widespread in RNA biology. The reg- 
ulatory expression platform (and the SD sequence) is part 
of the purine-rich region (Figure 1; orange). LI (green in 
Figure 1A) interacts with the major groove of the P2a/b 
helix in a triple helical structure that is characteristic of 
pseudoknot structures. The expression platform that is in 
the 3'-end of P2a/b follows the curvature of the helix and 
forms base-pairing with the 5'-end of P2a/b in the crystal 
structure of the metabolite-bound riboswitch. The nature 
of the unbound state of the SAM-II riboswitch at the 
atomistic detail is not known. However, the response of 
different regions of the SAM-II riboswitch to ligand 
binding has been extensively studied using chemical and 
inline probing experiments (10,39). Ligand-dependent sta- 
bilization of LI and the P2a/b helix were observed in the 
chemical probing experiments. On the other hand, the PI 
helix was predominately unchanged upon ligand binding. 

SAM binds in the groove between LI and the 5'-end of 
the P2 helix of the SAM-II riboswitch. The adenine group 
of SAM appears (Figure 1) to mimic a base on the 5'-end 



of P2 by intercalating and stacking between U21 and G22 
and interacting with its Hoogsteen base-pairing face to an 
unpaired U44 on the 3'-end of P2b, stabilizing the helix 
(Figures IB). U44 then interacts with U10 on LI, 
stabilizing the triple helix. The positively charged sulfur 
moiety of SAM interacts with the two carbonyl oxygen 
groups on U21 from P2b and Ull from LI, further 
stabilizing the triple helix between P2a/b and LI, in a 
site created by the triple base Ul 1»(U21-A45). This inter- 
action site seems to be one of the key components in the 
specificity of the SAM-II riboswitch for SAM, because 
SAH, which lacks the single positive charge, binds to 
SAM-II with an affinity reduced by at least three orders 
of magnitude (40). 

Atomic structures of the unbound state of riboswitches 
have been solved for only a handful of systems, and these 
structures include X-ray crystal structures of the glmS (33) 
and lysine (46,47) riboswitches and NMR structure of the 
preQj (48) riboswitch. NMR studies have also been 
carried out on the free and bound forms of the 
purine-sensing riboswitch (49-52). Most of the X-ray 
crystal structures of riboswitches that undergo large 
conformational changes have only been solved in their 
metabolite bound form. Therefore, questions do arise 
about the nature of their unbound states. For example, 
how do riboswitches switch conformations from the off to 
the on state, or vice versa? Direct observations at this level 
of detail are not always easy to achieve with current exper- 
imental techniques. In this regard, all-atom molecular 
dynamics (MD) simulation has proven invaluable as a 
complimentary technique to experiments in understanding 
biomolecular structures, dynamics and function (53-55). 
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Figure 1. The classic H-type pseudoknot structure of the SAM-II riboswitch (PDB id 2qwy). LI (green) is localized in the groove of the P2b helix, 
and the expression platform which follows the curvature of the P2a/b helix is part of the purine rich sequence shown in orange (A). SAM (cyan) is 
bound in the groove formed by LI and the 5'-end of P2b, and the adenine moiety of SAM mimics a base on the 5'-end of P2b by intercalating and 
stacking between U21 and G22 (B). The secondary structure of the SAM-II riboswitch is shown (C), with the base pairing interactions represented 
using commonly used notations. The other non-standard interactions observed in the X-ray structure that form the remaining tertiary interactions 
are not shown. The base-pairing interactions between LI, P2b and SAM in the metabolite bound X-ray crystal structure of the SAM-II 
riboswitch are shown in (B). The distances between N3 of U10 and 04 of U44 (black), 04 of Ull and 04 of U21 (red), and 04 of U12 and 
N6 of A46 (blue) are used to monitor the interaction of LI with the major groove of the helix. 
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This simulation method and related techniques have been 
used to study other types of riboswitches and have 
provided very useful information about ligand-induced 
conformational changes (56-59). We have therefore used 
MD simulations to probe the dynamics of the SAM-II 
riboswitch in the bound and unbound forms in order to 
understand the role of SAM in stabilizing the off state of 
the riboswitch and the nature of the expression platform 
in the bound and unbound forms. 



RESULTS AND DISCUSSIONS 

In order to study the atomistic basis for the switching 
mechanism between the on and off states of the SAM-II 
riboswitch, we have carried out two very long all-atom 
MD simulations (200 ns each) on the bound and 
unbound forms in explicit water and ions. The snapshots 
of the first 50 ns of each simulation were discarded in 
order to make sure that the systems were well equilibrated. 
The trajectories of the last 1 50 ns of each simulation were 
analyzed. 

Conformational switching of the pseudoknot structure 
of the SAM-II riboswitch 

There are three major hydrogen bonding and electrostatic 
interactions between SAM and the RNA (Figure IB) that 
are intricately involved in stabilizing the pseudoknot 
structure. The extended structure of SAM also spans 
three of the interactions made by LI and the major 
groove of P2b, as shown in Figure IB. We have monitored 
the three interactions formed by LI and the major groove 
of the P2b helix in the presence and absence of SAM in 
order to probe the role of SAM in preserving these inter- 
actions and the integrity of the SAM-II pseudoknot struc- 
ture. Two of the three interactions LI makes with the 
major groove of the P2b helix that are between N3 of 
U 10 and 04 of U44 and between 04 of Ull and 04 of 
U21 directly involve SAM and are well-formed in the 
complex (off state) throughout the simulation, as shown 
in Figure 2 (left column). The distance between 04 of Ull 
and 04 of U21 is used to monitor the second interaction, 
because the proximity of these two atoms is essential to 
the formation of the interaction site for the positively 
charged sulfur group in SAM. The third interaction 
between 04 of U12 and N6 of A46, also shown in 
Figure 2 (left column), is quite stable throughout the sim- 
ulation with some breathing motions taking place. The 
third interaction does not form any hydrogen bonding 
interactions with SAM. SAM locks the first two interac- 
tions in place and stabilizes the pseudoknot structure. The 
third interaction seems to naturally follow. 

We have studied these same interactions in the unbound 
form of the SAM-II riboswitch, also shown in Figure 2 
(right column). The interactions behave quite differently in 
the unbound form of the riboswitch. The interactions are 
not as stable as they are in the bound form. They form and 
break, alternating between different states. As a result, LI 
moves in and out of the major groove of the P2a/b helix. 
One of the states of the unbound form is close to the 
conformational state of the bound form, with some of 



the key interactions formed, and the other states can be 
viewed as a more relaxed and less compact ensemble of 
conformations. These results are consistent with chemical 
probing experiments of the SAM-II riboswitch (10). In the 
chemical probing experiment, the reactivity of a particular 
domain is directly correlated with its flexibility. It was 
shown that LI becomes less reactive upon binding SAM, 
suggesting that LI is flexible in the unbound form of the 
riboswitch and localized in the bound form. Distributions 
of the first two interactions that are made by LI with the 
major groove of P2b in the simulation of the unbound 
form of the riboswitch can be separated into a small pop- 
ulation that is close to the bound form and a more relaxed 
and open ensemble of conformations (Supplementary 
Figure SI). More contacts similar to that of the bound 
state are formed between U10 and U44 than between 
Ull and U21, and more contacts are formed 
between Ull and U21 than between U12 and A46 
(Supplementary Figure SI). Switching between the differ- 
ent conformational states of the riboswitch in the 
unbound form is estimated to be very fast and in the nano- 
second timescale (Supplementary Figure S2). Therefore, 
the limiting step of the on-off switch will most likely be 
the binding and dissociation of SAM that is expected to be 
much longer than the nanosecond time scale, based on 
binding studies of other riboswitches (60,61). 

These results suggest that the unbound form of the 
riboswitch periodically forms contacts similar to that of 
the bound form that are necessary to stabilize the 
ensemble of conformations of the bound form. The 
unbound form of the SAM-II riboswitch partially 
samples the conformational space of the bound form, 
and upon binding of SAM, the classic pseudoknot confor- 
mation is fully formed after 'zipping' up LI with the help 
of the additional RNA-ligand interactions. The simula- 
tion results therefore suggest that the recognition mecha- 
nism of the SAM-II riboswitch is due to a combination of 
conformational changes of the unbound form to confor- 
mations similar to that of the bound form and subsequent 
rearrangement of the structure, as exemplified by the 
conformational selection (62-64) and induced fit (65) 
models of ligand binding, respectively. Interestingly, this 
interplay between induced fit and conformational selec- 
tion has been previously observed in the guanine-sensing 
riboswitch (49). 

We have shown that the binding of SAM to the SAM-II 
riboswitch localizes LI in the major groove of the P2b 
helix, and our results provide atomistic description of 
earlier experimental observation (10) of the flexibility of 
LI in the unbound form. What functional role does LI 
have, since it behaves so differently in the bound and 
unbound forms of the SAM-II riboswitch and is a 
central feature in stabilizing the classic H-type 
pseudoknot? For example, how does the motion of LI 
correlate with other parts of the riboswitch? In order to 
answer these questions, we have calculated the 
cross-correlation coefficients of each residue with all of 
the other residues of the SAM-II riboswitch in the 
simulation of the unbound form, as shown in Figure 3. 
Figure 3 shows the normalized correlation matrix plot of 
both positive (red) and negative (blue) correlated motions 
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Figure 2. Distances between N3 of U10 and 04 of U44 (black), 04 of Ul 1 and 04 of U21 (red) and 04 of U12 and N6 of A46 (blue) in the bound 
form (left column) and unbound form (right column) of the SAM-II riboswitch during the simulations. The colors are consistent with the corre- 
sponding interactions shown in Figure IB. 
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Figure 3. Correlated motions of the SAM-II riboswitch calculated 
using the simulation of the unbound form. The color bar shows the 
relationship between the cross-correlation coefficients and the intensity 
of the colors. Red and blue represent strong positive and strong 
negative correlations, respectively. The residue numbers of the LI 
region of the SAM-II riboswitch that is in the major groove of the 
P2a/b helix are highlighted with a green line, and the residue 
numbers representing the purine-rich region of the expression 
platform are highlighted with an orange line. The positions of the 
two residues (U21 and G22) in the binding site, sandwiching the 
adenine moiety of SAM in the bound form are shown with arrows. 
The region showing correlated motions between LI and the expression 
platform is labeled as (A), and the region showing correlated motions 
between the binding site residues and the expression platform is labeled 
as (B). 



between each residue in the unbound form of the SAM-II 
riboswitch. Several positive correlations will naturally 
show up between different residues, if they are on the 
same domain or directly interacting with each other. For 
example, the diagonal shows up as positive correlated 
motions because it represents the correlation between 
each individual residue to itself. However, the interesting 
correlated motions are the ones that are long-ranged and 
do not involve direct interactions between the residues 
in question. For example, in Figure 3 we clearly see 
strong negative correlated motions between LI and the 
purine-rich sequence of the expression platform. We also 
see correlated motions between binding site residues, espe- 
cially residues U21 and G22, and the expression platform. 
These results therefore suggest and provide evidence that 
binding and dissociation of the ligand, coupled with the 
ensuing motions of LI in and out of the major groove of 
the P2b helix, are tied to the motions of the expression 
platform. The correlation of the motions between the 
binding site residues, LI, and the expression platform is 
interesting, because the mechanism of action of the 
riboswitch is to control the dynamics and conformational 
changes in the expression platform through binding and 
dissociation of the ligand. 

Ligand-induced curvature and base-pairing of the 
expression platform 

The expression platform is a key component of bacterial 
mRNA structure and directs the ribosome to the proper 
location in order to initiate protein synthesis. The ulti- 
mate goal of conformational switching in the SAM-II 
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riboswitch is to initiate or terminate translation due to 
induced conformational changes in the expression 
platform that expose or conceal the ribosomal recognition 
sequence. The AUG start codon, which is not part of the 
crystal structure, would normally be found downstream, 
immediately after the 3'-end of the riboswitch. As we have 
already seen above, the motions of the binding site, LI 
and the expression platform are correlated. A hypothesis 
based on occlusion of the recognition SD sequence in 
repressing translation was proposed as a result of the 
pseudoknot structure of the SAM-II riboswitch and 
chemical probing experiments (10). It was assumed that 
stabilization of the pseudoknot structure by SAM locks 
the expression platform in an ensemble of conformations 
that prevent the bases of the recognition sequence from 
ultimately interacting with the ribosome. One main obser- 
vation of the expression platform in the bound form (as 
shown in Figure 1) is that it is an integral part of the 
helical structure formed by the P2a/b helix, following 
the curvature of the helix, in the bound structure. Since 
the metabolite-bound form of the pseudoknot structure of 
SAM-II represses translation, it is obvious that this par- 
ticular conformation would not allow ribosomal assembly 
or translation to proceed. It is therefore assumed that a 
different ensemble of conformations of the expression 
platform would serve as a proper recognition conforma- 
tion or a proper conformation to initiate translation. We 
have therefore studied the curvature of the expression 
platform in the bound and unbound forms of the 
riboswitch, as shown in Figure 4. In order to calculate 
the effective curvature, we fitted the phosphorus atoms 
of six of the last seven phosphate groups to an effective 
circle. 

Figure 4 shows distributions of the effective curvature 
of the expression platform during the entire simulations 
with and without the presence of SAM. It can be seen that 




effective curvature,K (A 1 ) 

Figure 4. Probability distributions, pQc), of the effective curvature, k, of 
the expression platform in the bound (black) and unbound (red) forms 
of the of the SAM-II riboswitch. Ensembles of structures of the bound 
state (right inset) with a curved expression platform and unbound state 
(left inset) with a less curved expression platform are shown as insets. 
The structures represent a 5 ns segment of the simulation that matches 
the conformations at the peaks of the distributions. 



the expression platform of the bound form is more curved 
and has a narrower distribution than that of the unbound 
form of the riboswitch. The bound riboswitch exhibits a 
stable, tightly curved binding platform; however, the 
unbound riboswitch shows variation in curvature with 
an overall unwinding of the platform. The distribution 
of the effective curvature of the unbound riboswitch 
partly overlaps with that of the bound form, providing 
further evidence that the unbound form of the riboswitch 
partially samples the conformational space close to that of 
the bound form. The results suggest that in the absence of 
SAM, the platform uncurls and exposes the binding 
platform for possible ribosomal assembly or initiation of 
translation. As the binding platform straightens out, the 
bases of P2a/b become exposed and readily accessible, as 
shown in Figure 5A. 

Figure 5 shows the distribution of the interactions 
between the three most upstream bases of the purine-rich 
sequence (Figures 1 and 5B; orange) on the P2a/b helix 
that contains the expression platform on the 3'-end and 
the bases on the 5'-end of the same helix in the bound and 
unbound forms of the SAM-riboswitch. It can be seen that 
the interactions between these bases and their pairing 
partners on the opposite strand, as depicted in the X-ray 
crystal structure (Figure 5B), are very stable and well 
formed in the ligand bound form of the SAM-II 
riboswitch. However, in the unbound form of the 
riboswitch, these same interactions are loosely formed 
and the distributions are much broader than that of the 
bound form. Unlike the P2a/b helix, bases pairing in the 
PI helix are unaffected by metabolite binding (Figure 5C), 
also in agreement with chemical and inline probing 
experiments (10,39). Figure 5C shows distributions of 
the interactions between the bases of three base pairs 
around the center of the PI helix in the metabolite- 
bound and unbound forms of the SAM-II riboswitch. 
The distribution of the base pairing distances in the 
unbound form is almost identical to that of the bound 
form. 

These results therefore suggest pairing of the bases in 
the expression platform upon binding of the metabolite, 
concealing the bases from direct interaction with cellular 
components in the off state. On the other hand, the bases 
are hardly paired in the unbound form, and these loosely 
paired bases would be pivotal for the assembly of the 
ribosome. Alternatively, the expression platform of the 
metabolite-bound pseudoknot structure of the SAM-II 
riboswitch could still bind the ribosome, but blocks trans- 
lation until after dissociation of the metabolite and the 
subsequent conformational changes in the expression 
platform. Previously, Marzi et al. (66) reported a series 
of cryo-electron microscopy images suggestive of the fact 
that folded 5' UTR of mRNAs may still bind the 
ribosome, but blocks initiation of translation until 
the mRNA unfolds. The simulation results suggest that 
the changes in the curvature and base-pairing of the 
expression platform partly control the interactions of the 
riboswitch with the ribosome. 

Fully understanding how such conformational changes 
accompany ligand binding and how they direct initiation 
or termination of translation (or transcription, as in the 
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Figure 5. Probability distributions of distances, p(f), between the bases of three base-pairs in the purine-rich expression platform of the P2a/b helix. 
The distances are measured between N3 of U20 and Nl of A46, N6 of A19 and N3 of A47, and 02 of U18 and N6 of A48 (A) from the simulations 
of the bound (black) and unbound (red) forms of the SAM-II riboswitch. The base-pairing interactions are shown in the X-ray crystal structure of 
the SAM-II riboswitch (B). Probability distributions of distances, p(r), between the bases of three base-pairs around the center of the PI helix are also 
shown (C). The distances are measured between Nl of G3 and N3 of C29, N3 of C4 and Nl of G28, and Nl of G5 and N3 of C27 from the 
simulations of the metabolite bound (black) and unbound (red) forms of the SAM-II riboswitch. 



case of other riboswitches) would provide a valuable 
insight into bacterial gene regulation. The ability to 
modulate gene regulation in bacteria holds tremendous 
promise in designing new classes of antibiotics. The 
metabolite sensing riboswitches have so far only been dis- 
covered in prokaryotes, fungi and plants (13,67,68), there- 
fore, targeting this gene regulation mechanism as a 
therapeutic strategy could potentially have little or no 
side-effect to a human host. 



CONCLUSIONS 

Most of the X-ray crystal structures of riboswitches solved 
to date that undergo large conformational changes are in 
the metabolite-bound form. Therefore, there is limited 
structural information about their unbound states and 
the mechanism of the on-off switch. Our MD simulation 
studies have revealed the nature of the ligand-induced 
conformational changes of the SAM-II riboswitch at 
atomistic detail and agree very well with previous 
chemical and inline probing experiments. The classic 



H-type pseudoknot structure of the SAM-II riboswitch 
is only fully formed in the metabolite-bound structure. 
In the unbound structure, the SAM-II riboswitch alter- 
nates between a loosely formed pseudoknot structure 
and an ensemble of conformations that are close to that 
of the bound form. Ligand binding stabilizes LI in the 
major groove of the P2a/b helix, causing the helix to par- 
tially wrap around LI and induces curvature and stabilizes 
base-pairing in the expression platform. The loose and 
flexible nature of the unbound form of the SAM-II 
riboswitch would make it difficult to investigate the 
unbound form of riboswitches using X-ray crystal- 
lography. However, other experimental techniques, such 
as NMR, SAXS and chemical probing, could provide 
valuable insight into the nature of the unbound states of 
riboswitches. Similarly, theoretical and computational 
methods can complement experiments in further under- 
standing and providing another perspective of the 
switching mechanism of riboswitches in bacterial gene 
regulation at atomistic detail, as we have done in these 
studies. 
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Our simulations show that the rate of sampling the dif- 
ferent conformational states of the unbound form of the 
SAM-II riboswitch, including conformational states 
similar to that of the bound form, is on the nanosecond 
timescale. The results therefore suggest that the rate of 
dissociation of SAM will determine whether the mecha- 
nism of action of the SAM-II riboswitch is thermodynam- 
ically or kinetically driven. Since the lifetime of 
prokaryotic mRNAs is very short (minutes) (69), design- 
ing new drugs that can compete with SAM and bind to the 
SAM-II riboswitch with a very slow dissociation rate 
would therefore be a good strategy in shutting down the 
translational machinery of the SAM-II riboswitch. 

Computational methods 

MD simulations. The bound crystal structure of the 
SAM-II riboswitch (PDB id:2qwy; chain A) solved by 
Gilbert et al. (10) was the starting point for all simula- 
tions. All MD simulations were carried out using the 
Amber 10 suite of programs (70). In order to obtain the 
unbound form of the riboswitch, the structure was edited 
to remove the ligand. All of the crystallographic water 
molecules and ions were preserved. The modified nucleic 
acid parameters (parmbscO) (71) of the Cornell et al. (72) 
Amber force field were used. These parameters along with 
the structure of the unbound riboswitch were loaded into 
the LEAP module in Amber. The system was placed in an 
explicit TIP3P (73) truncated octahedron water box con- 
sisting of 1 1 257 water molecules, and 50 Na + counter 
ions were added to bring the system to neutrality. After 
minimizing the system, a 10 ns equilibration was carried 
out using the PMEMD module. The Langevin thermostat 
was used during the equilibration in order to bring the 
system from 100K to 300K using a collision frequency 
of 1.0 /ps. Following the 10 ns equilibration, a 200 ns 
long MD simulation was carried out. Particle mesh 
Ewald (74) method was used to treat long range 
electrostatics. A 2 fs time-step was used to solve 
Newton's Equations of Motion at a constant pressure of 
1 bar and a constant temperature of 300K using the NTP 
ensemble. All bonds involving hydrogen atoms were 
restrained using the SHAKE algorithm (75), and a 9 A 
cutoff was applied to all long ranged non-bonded 
interactions. 

Starting with the SAM ligand previously removed from 
the crystal structure of the riboswitch, the GaussView 
program (76) was used to add hydrogens to the ligand 
and the structure was quantum mechanically optimized 
using the Gaussian program at the HF/6-31-g(d) level of 
theory. Force field parameters were obtained using the 
Antechamber module in Amber, and partial charges 
were generated with the RESP (77) method from the 
electrostatic potential calculated using the Gaussian 
program, also at the HF/6-31g(d) level of theory. These 
force field parameters and charges were loaded into LEAP 
along with the default nucleic acid parameters in order to 
prepare the SAM-II riboswitch complex. The complex was 
also solvated in a truncated octahedron TIP3P water box 
consisting of 1 1 254 water molecules. The positive charge 
of the ligand resulted in the addition of only 49 Na + 



counter ions to the system in order to bring it to neutral- 
ity. Equilibration and dynamics were carried out from this 
point on under the same conditions and for the same 
length of time (200 ns) as for the free riboswitch. 

Cross-correlation coefficients. The residue-residue 
cross-correlation coefficients were calculated using the 
trajectories of the entire simulation of the unbound form 
of the SAM-II riboswitch. The cross-correlation coeffi- 
cient Cy = < Ar,Af7 >/(< Arj >< Arj >) 1/2 for the dis- 
placement of any two atoms i and j, where Ar, is the 
displacement from the average position of i. The 
cross-correlation coefficients, C,y, is therefore a 2D 
matrix with the value of each coefficient between +1.0 
and —1.0, where + 1.0 signifies perfect positive correlated 
motions, and —1.0 signifies perfect negative correlated 
motions. The cross-correlation coefficients were averaged 
over each residue, such that there was only one Cy per 
residue. The cross-correlation coefficients are then 
plotted using a color-coded 2D matrix from red (+1) to 
blue (-1). 
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